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Abstract 

In this paper difTerent types of compositions involving independent 
fractional Brownian motions Bjj, (t), t > 0, j = 1,2 are examined. 
The partial differential equations governing the distributions of Irit) = 
B}j,{\Bl^_{t)\), t > and Jpit) = B}j,{\Bl^{t)\^^"^), t > are derived 
by different methods and compared with those existing in the literature 
and with those related to B^ (\Bjj,^(t)\), t > 0. The process of iterated 
Brownian motion Ip{t), t > is examined in detail and its moments are 
calculated. Furthermore for Jp~^{t) = Bjj{\Bjj{. . . \B'Jj{t)\^/" . . .)\^^"), 
t > the following factorization is proved Jp^^{t) — Y[j=i (^)i ^ > 0- 
A series of compositions involving Cauchy processes and fractional Brow- 
nian motions are also studied and the corresponding non-homogeneous 
wave equations are derived. 

Keywords: Fractional Brownian motions, Cauchy processes, Modified Bessel 
functions, Iterated Fractional Brownian motion, Mellin transforms. Fractional 
equations. 

AMS: Primary 60J65, 60J60, 26A33. 

1 Introduction 

In the recent probabilistic literature there are some papers devoted to the inter- 
play between various forms of compositions of different processes and the partial 
differential equations governing their distributions. 

The best known example of composition of processes is the iterated Brownian 
motion (IBM) defined as 

Iit) = B'{\B'{t)\), t>0 (1.1) 
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where and are independent Brownian motions. The IBM has been in- 
troduced by Burdzy Q and its properties hke the iterated logarithm law, the 
fourth-order variation and many others have been analyzed in a series of papers 
dating back to the middle of Nighties (see Burdzy Q , Khoshnevisan and Lewis 
Q). The study of iterated Brownian motion has been stimulated by the analysis 
of diffusions in cracks (see for example, Chudnovsky and Kunin [S], DcBlassic 
0). It is well-known that the distribution of (|l.ip . say q{x,t), satisfies the 
fractional differential equation 

subject to the initial condition 

q{x,0)^S{x) (1.3) 
as well as the fourth-order p.d.e. (see DeBlassie Q, Orsingher and Zhao 16]) 

dq 1 d^q 1 



dt 23^4 2^2^ dx^ 



5{x), a; e R, t > 0. (1.4) 



Clearly -^^{x) = 5" must be understood in the sense that for every test func- 
tion < <j),5" >~ (/)"(0). This p.d.e. connection is vahd for any process 
X{\B{t)\), t > where X is a Markov process independent from the Brownian 
motion B (see AUouba and Zheng [l|, Baeumer et al. Q, Meerschaert et al. 
[loj . Baeumer and Meerschaert 0]). The time- fractional derivative appearing 
in (|1.2p must be understood in the sense 

d-'q,^^, 1 /•* ^q{x,s) 



and 



— — (x,t) = — / 7 — -; ds, m — l<v<m. 



d^q 

where m — 1 = \ . 

For the n-timcs iterated Brownian motion 

r{t) = B\\B\...{\B-+\t)\)...)\) (1.5) 
the distribution function solves the fractional equation 

xeM,t>0. (1.6) 

Extensions of these results to the d-dimensional case are treated in Nane (ll| . 
Orsingher and Beghin [l^ . The aim of this paper is to present different types of 
iterated processes constructed with independent fractional Brownian motions 
and to explore the partial differential equations governing the corresponding 
distributions. 

By considering fractional Brownian motions and Cauchy processes (suitably 
combined) wc obtain qualitatively different equations including a non-homogeneous 



2 



wave equation. In the probabilistic literature second-order hyperbolic equations 
emerge in the study of the telegraph process and of planar random motions with 
infinite directions. We show here that processes like C^(|C^(i)|) {C^ , j = 1,2 
are independent Cauchy processes) also are related to wave equations (see also 
Nane [12|). We consider here different types of combinations of processes in- 
volving the fractional Brownian motion and the Cauchy process. In particular, 
we show that for 

IcBAt)^Ci\BHm, i>0 (1.7) 

(where C is a Cauchy process independent from the fractional Brownian motion 
Bh) the probability law is a solution to the non- homogeneous heat equation 

d,_2m--^ .eM,t>0, (1.8) 



dt nx^^ 9x2 ' 

with H £ (0,1). In the special case H ^ ^ the p.d.e. (|1.8|) becomes the 
governing equation of the process 

IcBit) = C{\B{t)\), t>0. (1.9) 
For the vector process 

r^cit)^{B\\Cm,...,B-{\Cm), t>0 (1.10) 
we have that the joint law q ~ q{xi, . . . , a;„, i) solves the fourth-order equation 

= -^A2g-i-A5(xi)...5(x„), xeW\t>0. (1.11) 

This fact has been noted by Nane in fl^, [l^. It seems that it is difficult to 
obtain the equation corresponding to 

{B},^{\Cm,---,B^i^{\Cm), t>0 (1.12) 

where the independent fractional Brownian motions Bjj,{t) are involved. 
Finally we prove that the composition of two independent Cauchy processes 

Ic{t)=Ci{\C2m, t>0 (1-13) 
has a probability law satisfying the non-homogeneous wave equation 
d'^q 2 d-^q 



xeM, i>0 (1.14) 



which, in our view, is the most striking result of all these combinations of well- 
known processes. 

We also prove that the one-dimensional distribution of (|1.13p coincides with the 
multiplicative process 

Ci(|C2(<)|) ^Ci(y2t)C2(v^), i>0. (1.15) 

It is also true that the composition of Cauchy processes satisfies the following 
curious relationship 
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We remark that if the guiding process is a Cauchy process the related equation 
has a forcing function defined on the whole space-time domain, while in the 
case of a guiding process coinciding with Brownian motion, the forcing term is 
concentrated on the x line. 

The core of our paper concerns the iterated fractional Brownian motion whose 
general form is given by 

lF{t) ^ B]j^{\BlM), t>0 (1.17) 
where the processes B-'j^, (t), t > have covariance function 

E{B^^^is)B^^^it)] = 1 (|tp^^ + \s\'"^ - |t- sp^O > J = 1>2 (1.18) 
and < Hi,H2 < 1. 

For the iterated fractional Brownian motion 

init) ^ B\\Bj,m. i>0, i7e(0,l) (1.19) 

where is a standard Brownian motion independent from the fractional Brow- 
nian motion Bjj, we show that the equation governing its distribution q = q{x, t) 
is 

with q{x, 0) = S{x). 

The introduction of the fractional Brownian motion makes the equation (|1.20p 
slightly different from (|1.4I) in that it has non-constant coefficients. In the case 
where (|1.19p is replaced by the related process 

Lpit) = t'^BWBjj{t)\) t > 0, A' > 0, iJ e (0, 1) (1.21) 

the corresponding equation becomes a bit different and has the form 

(1.22) 

with an additional lower-order term with space-dependent coefficients. 
For the general iterated fractional Brownian motion 

lF{t) ^ B}j^{\Bj,^_{t)\), i >0, 0<i7i,i/2 < 1 (1.23) 

we have arrived at a substantially different equation which shows the funda- 
mental role of the process Bh^ ■ The equation pertaining to process (|1.23p reads 

(l + i/.i/.)t|+t^gf = i/,^i/.^{2x|+x^g}, xeR,t>0. (1.24) 

We note that equation (|1.24p has almost the same structure of (5.2) of Orsingher 
and Dc Gregorio [l^ which emerges in the study of random motions at finite 
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velocity in hyperbolic spaces. 

For the process (jl.23p we have the general expression of even-order moments 

which for Hi = H2 ~ 1/2 yields 

E{B\\B-mV'-^^^:^t''^' (1-26) 

which is the first formula of the remark 3.3 of Orsingher and Beghin 
Similarly for the n-times iterated fractional Brownian motion we have that 

(1.27) 

forn > 1, e N, ffo = 1 and Hj € (0, 1) for j = 1, 2 . . . , n . 
Another process obtained by composing independent fractional Brownian mo- 
tions is 

jUt) ^ Bj,{\Bjjit)\i), t>0, i/e(0,l) (1.28) 

which has proved to be much more tractable than (|1.17p . For (|1.28p it is even 
possible to obtain the explicit law in the following form 

f°° P~fl^ P~2t^ 1 /ItI\ 

(1.29) 

where K{){x) is the Bessel modified function defined as 
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Ko{x) = / -e^f^^'^ds, a; 7^ 0. (1.30) 
The p.d.e. governing the distribution p.29p is 

li'lM = -He--^{2^^^.^),M (1.31) 

= ^Hf"-^—{x — q{x,t)\, l-eM, t>0. 
We have also considered an extension of (jl.28p by taking into account the process 

Jl{t)^B\j{\Bl{\BUt)\^)\H t>0,He{0,l) (1.32) 

whose distribution is related to a fourth-order p.d.e. with non-constant coeffi- 
cients which reads 

do f 5^ (9'* 1 

J = i7t--|4— + 5.^ + .^^),, xeM,t>0. (1.33) 
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We note that for the extension of (|1.28|) . say Jp{t), the following equivalences 
in distribution hold 

,mt) = B},i\BU- ■ . is^' wi* ■ • oi^) =■ B\\B'i. . . \B"^+\t)f . . on, t > o. 

(1.34) 

Wc will also prove that for Jp~^{t), t > the factorization 

n 

r^-\t)=-Y[B'^it), t>0,He{0,l) (1.35) 

i=l 

in terms of independent fractional Brownian motions , * = 1,2 ... ,n holds. 
In view of decomposition (jl.35p one can also observe that 

Bh{\Bh{t)\^) = B\\Bj,{t)\')'d- B]^(t)Bl{t), t>Q,He{0,l]. (1.36) 

2 2 

We have considered some multidimensional compositions of Brownian motions. 
In Orsingher and Bcghin [l^, Nane [lH it is proved that the n-dimensional 
vector process 

{B\\Bm,...,B-i\Bm), t>0 (1.37) 

(where B^,...,B",B are independent Brownian motions) has a distribution 
function which is a solution to 

-dt^ = ^T.d^^ u{xi,...,x^,0) = l[S{x,), a;eM",i>0. (1.38) 

In AUouba and Zheng [l[ , DeBlassie it is shown that the law of (jl.37p satisfies 
the fourth-order equation 

S = ^fe£l - + :ri=±^,f[^i-rl .eR",i> 0.(1.39) 



\ J — 1 / J — 1 r—l 

For the n-dimcnsional process {B^{\BH{t)\), . . . , i?"(|i?H(i)|)) we arrive at the 
fourth-order p.d.c. 

g = i^=An^(-,) + ^^^A2g, .xeM",i>0 (1.40) 



which includes (|1.39p as a special case for = i. 

2 The iterated fractional Brownian motion 

For the iterated Brownian motion I{t) = B^ {\B'^ {t)\) , t > it has been shown 
by different authors that the distribution 



p{x,t) = 2l —=—=ds, xeR, t>0 (2.1) 
'2ns \/2iit 
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with initial condition p{x, 0) = S(x) is the solution to the fourth-order equation 
dp 1 d'^p 1 



dt 23 dx* 2V2^ dx^ 



d{x), xeR,t>0 (2.2) 



(see for example Allouba and Zheng [ij, DeBlassie 0, Nane [13). 
The distribution (j2.ip is also the solution to the time-fractional diffusion equa- 
tion 

subject to the initial condition q{x,0) = S{x) (see Orsingher and Beghin [l3| . 

M)- 

We now show below that the law 

pix,t) = 2 ^==ds, xeR, t>0, H,,H2e{0,l) i2A) 

Jo V27rs2-Hi V27rt2«2 

the iterated fractional Brownian motion Ipit) = Bj^ {\B'j^^{t)\), i > is a 
solution to a second-order differential equation with non-constant coefficients. 
We first need the following auxiliary results concerning the density of fractional 
Brownian motion. 

Theorem 2.1 The fractional Brownian motion Bnit), t > 0, H £ {0, 1) has a 

probability density pH{x,t), x & t > Q which solves the heat equation 

with initial condition ph{x,0) — S{x). 

Proof: The distribution of Bnit). t > can be written as inverse Fourier 
transform 



PHi^,t) = ^ 1 e-^^^e-^^'''"d(3. (2.6) 



A simple calculation proves that (|2.5p holds. For H ^ ^ equation (|2.5p coincides 
with the classical heat-equation. 



Remark 2.1 For the Gaussian law 

q{x, t) = 4==' x£R,t>{) (2.7) 
^J2iig[t) 

with g £ , g : (0,oo) i— >■ (0, oo), we can easily accertain that the governing 
equation is 

d do 1 d'^ 

di'=i2dx^'^ x-eR,t>0. (2.8) 
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We prove below that the processes = B^{\B'jf{t)\) and its shghtly 

extended version Lnit) = B^{\Bjj{t)\) have a distribution satisfying a fourth- 
order p.d.e. somewhat similar to (|2.2p . We analyze in the next theorem the law 
of Lnit) and extract the related equation for I nit) by assuming that K = 0. 

Theorem 2.2 The law of the process 

Lpit) =t^BWBjf{t)\), t>0,K>0 (2.9) 
is solution of the fourth-order equation 

(210) 

subject to the initial condition q{x,0) ~ S{x). 

Proof: Wc start by writing down the distribution oi Lp{t), t > as 



^) = 2 / -_=-===ds, xeR,t>0,K>0. (2.11) 
Jo V27r.st2«; V27rt2« 

By taking the time derivative of (j2.1ip we have that 

dp f°° d ( e^i^ \ e^i^H- f°° e^"^ d I e^^ , , 



dt Jo dt \ V27rst2^ / \/2^rt2ii- V27rsi2if 5f \ 

(2.12) 




2i7t 



V2TTst^^ ds^ \ 



In the last step we used equations (|2.5p and p.Sp . The first integral in (|2.12|) 
can be developed as 




t2K g 2t2H K d , . 

zds) = — -— [xp) . 



V2^^^ V2^^ j t dx 
By integrating by parts the second integral in (j2.12D we have that 



2i/t 



2H-1 



y/2sTTt^i< ds I V2^rt2H- 



Jo 5s I V27rsi2^ ) ds \ 
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-2Ht 



2H- 



e 2t 



2Ht 



CO Q2 



s=0 



zds. 



By assuming cj) = st^^ and in view of (|2.8p we can write that 



d e 23*2 



All this leads to 

K d 



2K Ql ^ 

and 



(9p 
9f 



H 



-1 5^ 



e 23*2 



e 2t 



9a:;2 \ V27rsi2A I V2^rt27i- 







V27rst2^f y2^rf2H 



ds 



K d , , m^+2^-1 d2 

t OX 



2tt dx' 



dx'^' 



This concludes our proof. 



Remark 2.2 For K = 0, equation (|2.1ip simplifies and becomes 



t 



dq H^^^d^Q Ht" d2 



'dt 22 9a;4 dx2 
Clearly for H = 1/2 we reobtain DeBlassie's result 
dq 1 d'^q 1 ^2 



(5(a:), x e K, i > G (0,1). (2.13) 



dt 2^dx^ 2x/2^dx^ 



5{x), X e M, <> 0. 



(2.14) 



It is useful to have a look at table 1, for comparing the governing equations 
examined so far. 

It is easy to check that the vector process {B\j {t), . . . (<)), i > has a law 
which satisfies the heat equation 



dt 2^ ^ -^-2 



-id'^q 
9^2 



(2.15) 



For the n-dimensional vector process {B^{\BH{t)\), . . . B'^{\BH{t)\)) the distri- 
bution function can be written as 



q{xi,...,Xn,t) = 2 ds e »2.,=i".^j ^ ^^=^"^ ^==dai . . . dar, 
Jo JR" V27rt2H 
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Process 


Governing 


Iterated 


Governing 




equation 


process 


equation 


B{t) 


dq 1 d'^q 
dt ~ 2 dx^ 


Bmm 


99 _ 1 1 1 X 
dt ~ 2^ dx-i ^ ^jMdx-'" 

d^^^q 1 d^q 
dt^/^ ~ 23/2 


Bnit) 


dt dx^ 




+ dx-^°y^) 



Table 1 



(2.16) 



By taking the time derivative and considering (|2.5p we have that 



99 

di -Jo (2^)» 



. . dan 



(27r)'^ 



_^ d I e 2t2H 



dai . . . dan 



s=0 



-Ht 



2H-1 



Ht 



2H~1 



ds 

(2^ 
ds 



d 



e 2t^ 



dcki . . . da^ 



Ht 



2H~1 



(27r) 
1 



9 



e 2t^ 



as 1 



dai . . . dan 



-Ht 



2H-1 



(27r)" 
°° ds 



(27r)^ 



E 



V2^ rfa;^ 



22 



" ^2 \ / " ^2 



Ht 



2H-1 



2H 



\f2wt 



where A — y'"_i 



A(5(a;i)---(5(x„) + 



2ff-l 



22 



The iterated fractional Brownian motion l\(t) ~ Bjj_^ ( | Bj^^ (t ) | ) brings about 
substantiaUy different p.d.e.'s and we start our analysis from the first-order 
equation of the next theorem. For a duscussion on the uniqueness of Cauchy 
problems related to the above equations see Baeumer et al. 



10 



Theorem 2.3 The distrihution of the iterated fractional Brownian motions 
Ip{t) = B\j^ {\Bl^^{t)\) is a solution of the first-order p.d.e. 

dv d 

t-^ = -HiH2g^ixp), xeR,t>0 (2.17) 
where Hi, H2 e (0, 1). 
Proof: The time derivative of 



Jo V27r.s2^i V27rt2^f2 



writes 



p(x t) = 2 ^^=^= lis 

' Jo V27rs2Hi dt V27rt2ff2 



00 



^2 ^2 

e 2a2Hl ^ g 2*2^*2 



V27rs2^^i (9s V27rt2ff2 .=0 V27rs2ffi 9s V2^d^ 



d e 2b^"i d e 2t^"2 

zds 



io as V27rs2«i as V2^^ ^ 
-2H2t'-"^-' / i/is2«-i^4^_ f-^) ' ds 

5x2 ^/2^^ V ^2^2 J ,/2^ 

^2 ^2 

2HiH2t~^ / s^^^ ^ , ^ ^ zds 

io 5x2 ^27rs2ffi V2^rt2i?7 



= -2HiH2t-^^ X / ^^=ds = -HiH2t-^-^ (xp) 

dxy Jo V2^^ V2^- j dx^ 

and this concludes the proof of theorem. 



We study in the next theorem the ensemble of solutions to (|2.17l) 
Theorem 2.4 The first order partial differential equation 
t du du 

+ x— = -u (2.19) 



H1H2 dt dx 
has a general solution of the form 

u{x,t)='^f(-^y xeK-{0}, t>0 (2.20) 

with Hi, H2e (0, 1) and f e C^{R). 
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Proof: The auxiliary equations pertaining to (|2.19p are 

H,H,^^^^^^. (2.21) 
t X u 

The first couple of members of (|2.2ip 
t X 

has solution 

while equating the second and third term of (|2.2ip we get 

log x = — log u + constant 
and thus 

w = -C2 = -J (ci) = -/ 

XX X 



tHiH2 



Remark 2.3 We note that the functions of the form 

\ / X \ 1 / X 

u{x,t) = -f [jhTh;) = jHTH^a [jhTh; 

with f,g G Ci(K) are solutions to ((2Tf| or ([2T9)) . 

The law of the iterated fractional Brownian motion belongs to the class of 
functions (|2.20p because it can be written as 



or 



1 X e 2t^«i^2,2i^i _ 1 / X 



Jl ,2 

Furthermore if we choose /(z) = ^'^^j^ or g(z) = ^/=- we see that the Gaus- 
sian kernel of Brownian motion belongs to the class of solutions to (|2.17p with 
H1H2 = ^■ 

Theorem 2.5 The density of the iterated fractional Brownian motion Ip{t) = 
B},^i\Bl^m for HuH^e (0,1) 

^2 s2 

q{x,t)^2 ^==ds, .tgR, t>0, i7i,i/2e (0,1) (2.22) 

is a solution to the following second- order p.d.e. 

{l+HMtf^+t'^ = HfHl |2x^ , xeR,t>0, H,,H,e (0, 1). 

(2.23) 

We give now two different proofs of this result. 
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First proof: We start by taking tiie time derivative of (|2.22p and by taking 
into account (|2.5p we have that 



dp 

'di~~ Jo V27rs2-ffi 



^ d e 



dt \yV2^^ 



ds 



= 2H2t 



2H2-1 



ds 



V27rs2ffi ds^ \ V27rt2ff2 

die 2s2Hi \ d I e 2t^^2 



ds \ V27rs2ffi I as I V27rt2ff2 



ds. 



(2.24) 



In the last step an integration by parts has been carried out. 

The critical point emerges on earring out the second integration by parts which 

yields 



at 



^52 I V27rs2ffi / ^/2^^-' 



■ds 



(2.25) 



as 1 V27rs2«i / V27rt2ff2 
By applying once again (j2.5p we have that 



a7 = 2^^^ 



e 2s'^^i I e 2t^"2 
as2 I V27rs2ffi J Tl^P^ 



ids 



(2.26) 



^2 



c'a;^ \ V27rs2«i I V27r<2ff2 



s=0 



From (j2.26p it is clear that for Hi > \/2 the second term disappears for s — >■ 0+ 
while for < i?i < 1 /2 we get a Dirac function as in all previous cases examined 
above. 

It is convenient to rewrite (j2.26p in the following manner 



dp 
'dt 



^ = 2H-,f"^-^ 



00 /^2 



52 / 2,2Hi \ 2,2if2 2i?li72 R 1 . ^ 

l =1 ds+ ' ' t^^'4(oi/2l 

9s2 \ V27rs2ffi ) V2^rf27f7 ^ ("'i/^jv ;^^2 

(2.27) 



where the first integral can be developed as follows 

^2 



_5 

ds 



2i/2i'^^-^ / ^ <^ His'' 



dx^ \ V2^^^ I V27rt2Jf2 



ds (2.28) 
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and can be conveniently rewritten as 



Jo 

We note at this point that for Hi ~ 1/2 the first integral in the last member of 
(|2.28p disappears and immediately we obtain equation p.l3p . For Hi ^ 1/2 an 
additional analysis is necessary and a qualitatively different result is obtained. 
The time derivative (|2.27p multiplyed by t^~^^ becomes 



t^-"^^=2HiH, 
ot 



„2Hi-2 



_ ( e'A 



a e 



e 2t2«2 iJiiJo 

2, '"+^T71#^<«*("''S^*(^>' 




To get rid of the s ^ appearing in the above integral we derive both members 
w.r. to time t and realise that 



a_ 

m 



1-H2 



dp 
di 



(1 - Hi) 



~H-\ 



, HiHi 
't^^+i Jo 

dx^ \'" V27rs2«i 



e 2s2«i 



a 



aa; I V27r.s2ffi 



e 2=2^1 



e 2t2"2 
V27rt2ff2 ' 



H1H2 

t^2+l 
H1H2 

H1H2 
tH2 



dx 



'j\{l + Hi)—{xp)+Hi 
^^\i^ + Hi)§^{xp)+Hi 



a2 

aa;2 






a 






dx 


} 


2x 


dp , 2 
ox 


a^p 


ax2 



In light of theorem (|2.5p . by writing ~HiH2^{xp) = and after some cal- 
culations, equation (|2.23l) emerges. 
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Second proof: By taking the time derivative of (|2.17p multiplyed by t we 
have that 



d f dp 
*dt Vdt 



dp 2^ 

d f dp 



-H^H2-^{xp) 



HfHi\li.p)^.^i.p) 



= HlHl^[xp) + HlHl 



2x 



d_ 

dx 



= -HiH2t 



dp 



HlHl 



dx'^ 



2.t|^+x2 — 
ox ox^ 



,d^p 



In the last step we used again (j2.17[) . 



Remark 2.4 The equation (j2.23p displays the structure of the telegraph equa- 
tion with non-constant coefficients emerging in the case of finite-velocity one- 
dimensional motions in a non- homogeneous medium (see Ratanov [l3|)- An 
equation similar to (j2.23p emerges in the study of hyperbolic random motions 
of finite velocity in hyperbolic spaces (Orsingher and De Gregorio In view 

of remark 12.31 and also by direct calculations it can be proved that the distribu- 
tion densities of fractional Brownian motion, Brownian motion are solutions to 
the wave- type equation (j2.23p . 

For the n-timcs fractional Brownian motions we can easily evaluate the even- 
order moments. 

Theorem 2.6 For the n-times iterated fractional Brownian motion 

mt)^Bj,^{\Bj,^{...\B-+lit)\...)\) 
the explicit form of the moments reads 

2"+ii2fen"i^ff. " r(2fcn';=o^j) 



Emm 



2k 



2'=^"-n;.„H, n r(fcn;'=o^.) 



fc e N, 77. > 1 (2.29) 



with Hq = 1 and Hj G (0, 1) for j = 1,2 . . . ,n. 



Proof: We proceed by recurrence by first evaluating 



^^^^^ ' V2^^V2^^ 
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Since 



P 23^^! ok+i 



r(2fc' 

m 



we obtain that 

_ t2fcgig2 r(2fcgi) 

^ 2/=^! T{kHi) 

and thus 

P r ji c+u2fe ^ 22-fc(i+«i)f2fcHiif2 r(2fc)r(2fcgi) 

r(/c)r(fci/i) ' 

If we assume that result (|2.29p holds for n — 1 in view of (|2.30p we immediately 
have the claimed formula. 



Remark 2.5 The variance of the iterated fractional Brownian motion 

Var {B},^{\Bj,M)} = 2^-"^^-^^"^"^ (2.31) 

is an increasing function of the time t which grows more rapidly than the vari- 
ance of IBM 

Var{B\\B\t)\)]^^ (2.32) 
for all Hi,H2 such that 
H,H, > 1. 

For H1H2 > 1/2, Var {Bjj {\B'jj^{t)\)^ increases more rapidly than the variance 
of standard Brownian motion. 



3 The iterated fractional Brownian motion and 
modified Bessel functions 

A process related to the iterated fractional Brownian motion (i) , t > is here 
analysed. Unlike the n-times iterated Brownian motion 

I"{t)=B^(\B^{...\B''+\t)\...)\), t>0 (3.1) 
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for which the probabiUty distribution is expressed as a 71-fold integral, for the 
process 

JUt)^Bj,{\Bj,{...\B-+\t)\i...)\i), t>0 (3.2) 

the probability density can be explicitly expressed in terms of modified Bessel 
functions. In particular, for n = 1, 2, 3 we have the following expressions 



Pr { 4 (t) e dx} = dx-^Ko (^^^ 



(3.3) 



Pr {Jl{t) e dx] = dx -^__Xo (^^ ) ds (3.4) 



Pr{jm^dx)^dx [ ^K.(\^)^K.(\-^)ds. 



(3.5) 



In general, for n ~ 2m we have that 
Pr{jj^''{t) edx]^ 



dx I ds,... I ds,^—Ko ( M ) . . . ^— /vo ( ) -^Ko 



(3.6) 



and for n = 2m + 1 



Pr {Jl'^+^t) e dx} =dx ^=Pr{jp\t) edsi} (3.7) 

Jo v27rsj 

= / e C?^} , ^ dSr.. (3.8) 

Jo v27rt^^ 

We shall call J]?(t) the weighted iterated fractional Brownian motion. 
The function K^{x) can be presented in several alternative forms as 

Ko{x)= s'^e^f^^'^'ds, |arga;|<-- (3.9) 
Jo 4 

(see p. 119 Lebedev 0) or, for ^{v} > -i 

X.(.).(i^)^^^ re-^.^-^(l + A)^-d., larg.K. (3.10) 



2) "'O 

(see p. 140 Lebedev 0) or, for 7^ 0, u = ±1,±2, . . . 

^ , arga;|<7r (3.11) 

2 sm :/7r 

(see p. 108 Lebedev [§]) where 

/.(-) = E(2) fc!r(fc + . + l) ' \-\<oo,\-rgx\<n (3.12) 
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is the Bessel Modified function of the first kind (see p. 108 Lebedev Q). 
For the process 

jUt) ^ B},{\B%it)\i), t>0 (3.13) 
the distribution function, in view of formula p.9p . becomes 

(314) 

where the change of variable s = \/2zt^ must be introduced. 
In the spirit of the previous sections we give the partial differential equations 
governing the distributions p.3p and p.4p . We now state our first result for the 
process Jp{t), t > 0. 

Theorem 3.1 For the process 

Jhit) = B},{\Bj,{t)\^), t>0,i/e(0,l) (3.15) 

obtained by composing two independent fractional Brownian motions B\j and 
B^, the distribution iS.14^ for x ^ 0, solves the p.d.e. 

for all X ^0, t > and H e (0, 1). 

First proof: We derive equation (j3.16|) by two different methods. The first 
one repeatedly uses the relationship of theorem (j2.ip while the second one is 
based on the properties of the modified Bessel functions Ky{z). 
We start by calculating the time derivatives of (|3.14p as follows: 



= 9 



V27rs2 dt 



Ph[s, t)as 



s=0 



d / ert^ \ d , , , 

r ^ r- ^ -^PH[S,t)ds 

ds \ V27rs2 / ds 



Since 



s 

—Ph{x,s) = --^pH{s,t) (3.17) 

it is easy to see that both functions ^^^=f and -^pnis, t) tend to zero as s — >■ oo. 
For s — )■ clearly -^pnis, i) — while the other term converges to zero because 
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e — 7> as s for X 7^ 0. 
This means that 



e 2s 



= -2Ht 



ds \ V27rs2 / ds 
2H-1 J d I e-i 



PH{s,t)ds 



ds \ ^/2^ 
= [ by (EH)] 



PH{s,t) 



s=0 



oo Q2 



e 2s-! 



9s'^ \ V27rs2 



PHis,t)ds 



= 2Ht 



2H-1 



Ht 



2H~1 



9^ I 727^2 



/ e 2l^ 



a; e 23 



PH{s,t)ds I . 
PH{s,t)ds 



Hf-^~'lS^p\x,t)-^[xpHx,t)] 



In the previous steps it is important to consider that 



dx^ ' dx^ 
92 93 1 , 



d f 1 



-,a^_ia2,2^^\ = J- 



'd's\^/2^j dsy2nj_^~' 27r , 

where the inversion of derivative and integral is justified since 

122 



possesses finite integral. For the same reason the limit w.r. to s can be brought 
inside the integral and 



lim 



>o+ ds \ V27rs2 



PH{s,t)=0. 



Second proof: This proof is based on the following properties of the modified 
Bessel function (see p. 110 Lebedev [§]) 



and 



az z 



K,+i{z)^K,_i{z) + 2-K,{z). 



(3.18) 



(3.19) 
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We perform our calculations on the representation p.l4p of the distribution 
of the iterated fractional Brownian motion. We need to evaluate the partial 
derivatives appearing in equation p.l6p for x G R"*", for the derivative w.r. to 
time we have that 









{w)\ 





H 



Hx 



(3.20) 



by dSH] 

Ht 



H 



We now pass to the partial derivatives with respect to space. Wc have succes- 
sively that 



d_ 

dx 

&^ 

dx'^ 



(^) 
(J) 



by (ixra 

1 d 



(3.21) 
(3.22) 



Ki{z) 



1 



H 



-Kn ^ A'l 



X \ 



Sonic more calculations are necessary for the third order partial derivative which 
reads 



dx^ 



1 



Ko ( ^ 



^{^^^{iH) + '-^K,{^) (3.23) 



x \ 



— A'l _ - if, _ 



1 



by 



1 



1 



-A'o ( ^ ) + ( 2— + — I K, ( — 



Therefore, by combining the previous results one obtains 



8P_ 

dx^ 



1 



Ko ( JH 
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and thus 



+ (^^)] = JtT {-^0 ( 

which coincides with the time derivative p.20p . 
This concludes the second proof. 



)} 



We now examine the twice iterated fractional Brownian motion and show 
that its probability distribution solves a fourth-order p.d.e. 

Theorem 3.2 The process Jl{t) = B]j{\B'jj{\B^u{t)\H)\H) has a probability 
density 



V27rs2 Jo V27rz2 V2^d^ 
e 27^ 1 



zdzds 



y2/KS^ 



(^) 



(3.24) 



p\s,t)ds, X eR, t> 0, H e {0,1) 



satisfying the following fourth-order p.d.e. 
/or x 7^ 0, i > 0. 



(3.25) 



Proof: We start our proof by writing down the time derivative w.r. to time t 
of (|3.24p and making use of the result p.l6p we get that 



(3.26) 



-2Ht 



2H-1 



I— s — 

ds'^ ds-^ 



p^{s, t)ds. 



By applying a couple of integrations by parts we have that 



s ^ <^ 1/ \j I o I e 2s 
V2^d^^^''*' Jo 71^ 



p^{s, t)ds 



(3.27) 



while by applying three integrations by parts we obtain 

g3 



/ 

Jo 



6 2s'^ 



e 23 



9s3 1 V27rs2 



p\s,t)ds (3.28) 
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In light of p.27p and (|3.28p the time derivative p.26[) becomes 
d 



92 



92 



ds-^ \ V2^ / 



d / e: 



ds \ V27rs2 



■ p^{s, t)ds 



-2Ht 



2H-1 



2Ht 



2H-1 



92 

952 \ V27rs2 j ds 
92 / e- 



d I e ^ 



ds \ 7271 



g3 



+ s 



9s2 I V27rs2 / 9s3 \ V27rs2 
For simphcity, we now write 

g = qix,s) 



9s2 1 V27rs2 
p^(s,i)ds. (3.29) 



V27rs2 

By repeated use of (|2.8p we have that 



92 



9 



i9s2 9s 



92:2 



92 



r9 + s 



9x2 



and 



9^9 9 ^2 
9s3 9s 



9^2 9a;'' 



9x^ 9x^ 



Hence, by combining the previous results we can rewrite the time derivative 
(P:^ as follows 



d_ 

dt 



p2 (x, t) = 2i/i2H-i J J —qix, s) + 4s2 ^9(2;, s) + s) (s, t)ds 

(3.30) 



In order to handle the operator 



r_9^ 



+ 4s2 



94 



96 



\ 9a;2 ' 9a;^ ' " 9a;^ 

it will be useful to note that 

9 X 
—q(x,s) = -—q{x,s) 



qix,s) 



(3.31) 



and 



92 



1 x' 



-Q^q{x,s) =(-— + —) g(x,s), 
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this implies that the operator p.3ip becomes 



dx"^ 

dx"^ 
92 



9 + 4s' 



53 
dx^ 



X 



Q3 

4 

dx'^ dx^ ^ 



dx^ 



dx^ 
94 



dx^ 



[x\) 



q3 



By some calculations we have that 



93 
9a:3 



(x q) 



and 



54 



92 33 



5' 2d 

——^ 2xq + X —q 



dx^ 



9^2 



9^2 



q + X 



q + X 



dx 

93 
9x3 

93 
9x3 



9x2 



2x— g 

ox 



9x2 



9+ ^ 



9_ 

9x 



_9_ 

9x 



n 5' 2 5^ ■ 

ox^ 9x'^ 



12— 

9x2 ^ 



93 2 94 

ox"^ ax"* 



by inserting p.33p and p.34p into p.32p we have that 

92 2 94 4 96 92 93 

9x^ 9x* 9x° 9x^ 9x'' 

and this leads to the claimed result. 



9x2 



9x2 



(3.32) 



(3.33) 



(3.34) 



Let's consider now the process 
Jm = Bh^{\Bl^{\ . . . \B-tUt)\'^ • ■ • i > 0, iJ, e (0, 1), J = 1, 2, 



We have already seen a particular case in p.34p and (|1.36p where Hi = H2 = 
. . . = Hn+i = H . We show that the following equivalence in distribution holds 



Jm = B]j^{\Bl^{\...\Blll^{t)\^^ 



|*)|*)^^-Si(|i32(...|B«+;^(t)|^ 



(3.35) 
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It must be also noted that the unique relevant Hurst parameter H in Jp{t) 
is that in B^^_^{t). For this reason and for the sake of simplicity we assume 
throughout the paper that Hj = H for j = 1,2, . . . ,n + 1. 

We can state the following result on the factorization of the iterated fractional 
Brownian motion. 

Theorem 3.3 For the process 

Jr'it) = Bhi\BU- ■ ■ \Bm\^ ■ ■ Ol*), t>0,HeiO, 1) (3.36) 

with B^, j — 1,2, ... ,n independent fractional Brownian motions, the following 
equalities in distribution hold 

n 

)!')=■ n^fW' ^>0' He (0,1) (3.37) 

1=1 

where (t), t > 0, j = 1,2, ... ,n are independent standard Brownian motions. 



Proof: The proof of the first equality in p.37|) is a trivial matter. Indeed the 

equivalence in distribution follows from the property Bnit) B(t^^) for each 
t. For the proof of the second equality in p.37|) we proceed by evaluating the 
Mellin transform of the density of \Jp~^{t)\, i > as 



Mn-i{a,t) 



22 F I 



dsi . 



e 1 e 



e 



i-Sn-i- 



dx 



The Fourier transform of the distribution of ^ (t) can be written as 



J-„_i(/3,t) = / e'^''pn-i{x,t)dx 



E 



ki 



fe=0 



E 

fe=0 



m 



2k 



{2k)\ 



X^Pn-lix, t)dx 
X^^Pn-l(x, t)dx 

2^+^r(fc + i) 

\/2^ 



,2Hk 



We now evaluate the characteristic function of 

n 



which becomes 

i/3n;=iB'„(*) 

he 



fe=o ■ I j=l 



(3.38) 
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fc=0 j=l 



fc=0 ^ ' j = l 



2 A: 



^ (2fc)! 

fc=0 ^ 1 = 1 



2'=+^r(fc + i) 



/2^ 



E 

fe=0 



2fc 



(2fc)! 



2^+^r(fc + i; 

71^ 



,2/ffc 



and this coincides with p.38p . This confirms the equivalences in distribution of 



Remark 3.1 The result of the previous theorem permits us to write down the 
following special factorization 



BU\Bl{t)\^)'^- BUt)Bl{t)'^- B\\Bl{t)\'), t>0, ife(0,l). 



4 The iterated Cauchy process 

The Cauchy process has been investigated from many viewpoints and some au- 
thors have examined the structure of the jumps, some others have concentrated 
their analysis on the behaviour of sample paths and others have obtained a 
number of its distributional properties. We here consider the composition of 
independent Cauchy processes and its connection with wave equations. 
The Cauchy process C{t), t > has a law 



which satisfies the Laplace equation 
a2 d^T) 

+ xeR,t>0 (4.1) 

and also the space-fractional equation 

While (|4.ip can be cheeked straightforwardly some additional care is needed for 
(113). 

The first order spatial derivative appearing in (|4.2[) is defined as 
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We now evaluate the Fourier transform of (|4.3|) . 



TT Us {X - S) J„oo (S - X) 

TT \Jr 



It is "well-known that 

2 r-sin^^^^M 

Jo y p 

and thus the Fourier transform of the fraetional derivative (14.31) becomes 



e'^^l^l^) {x)dx^\f3\J^^f'^f{x)dx 



Since 



s-"^ ' ,. dx = e-*l^l (4.4) 
equation (j4.2p immediately emerges. 

The iterated Cauehy process defined as 

icit) = c'i\c'm, t>o 

has a probability law represented by the function 

s t 



2 /"^ 



ds X e M, t > 0. (4.5) 



+ a;^ + 
We prove that 

Theorem 4.1 T/ie function q = q{x,t) has the following explicit form 
2t t 

,(x,t) = -,^^^log- xeR,t>0 (4.6) 

and satisfies the non-homogeneous wave equation 
d'^q d^q 2 



dt"^ dx^ 'kHx'^ 



X e M, t > 0. (4.7) 
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Proof: The integral can be easily evaluated by observing that 
t st ( 1 1 



+ t'^ + 



+ 



The proof that q satisfies the equation (j4.7|) can be carried out in two alternative 
ways. 

By taking the second-order time derivative of (|4.5p and by applying (|4.ip we 
have that 



•7r2 , 



92 



t 



ds 



2 

~2 



°° d / s \ d f t 



ds 



ds 

2_d_ /_ 

TT^ ds V 



00 



ds 



32 2 



TT^X^t dx'^ TT^ 

2 

TT^X^t dx"^ 



92 



(9s2 V 5^ + a:;^ / 



rds 



-ds 



The second proof that (j4.6p satisfies (|4.7|) requires some additional attention 
and is based on the explicit calculations 



d\ _ 2H{e + 3x2) ^ 



log 



+ 3i2) 



^2(^2 _ 2.2)3 '=2; 7r2t(t2 _ a;2)2 



92g 22t(t2+32;2) f 

- log- 



(?a;2 7r2(i2 _ a;2)3 
which readily give that 

9^ _ 9^ 2_ 

at2 dx"^ TT^X^t 



7r2x2(t2 _ ^2)2 



Remark 4.1 The p.d.e. (|4.7p is part of a more general result of [12, Theorem 
2.1]. 

The iterated Cauchy process has an equivalent representation in terms of 
products of independent Cauchy processes as shown in the next theorem. 

Theorem 4.2 The following identity in distribution holds for all t > 



CWC^{t)\) ^C\V2i)C^{V2i). 



(4.8) 
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Proof: It suffices to show that for w > 













TT^ Jo Jo 



2t V2t 



and therefore the density p — p{w, t) of the random variable (vSi)^ (v2i) 
becomes 



2^t f 



X 



22t 



log- 



2t + x2 te2 + 2^2 

2t + a;2 oo 2t 



1 !^ 

7r2(u;2 _ ^2) ^2.2 2u,2 7r2(w2 - t2) °g i • 
Similar calculations must be done for w < Q. 



Composing more than two Cauchy processes make the research of their dis- 
tribution and of the governing equation a puzzling question. 

Two different compositions of the Brownian motion and Cauchy process are 
now analysed. We first obtain the p.d.e. (|1.8p which is the governing equation 
of the process IcBh ~ t > {). The transition law of the process 

Icbh is given by 



p{x,t)^2 



e 2t^ 



-ds. 



lo 7r(s2 + a;2) y^2^ 
By performing the time derivative of (|4.9p we obtain that 



dp 



d 



e 2t^ 



2m 



7r(s2 + x2) dt \ V2^ 
92 



ds 



2H-1 



7r(s2 4-a;2)as2 \ ^2^^ 



ds 



-2m 



2H-1 



d 



d 



as V7r(.s2 +a;2); as \ V2^^ 



2i7f 



ds V7r(s2 + a;2)y ^2^;^ 

oo ^2 



as2 V7r(s2 +a;2)y ^72^^ 



zds. 



(4.9) 
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It must be observed that the first term of the sum becomes 

1 1 2Ht"-^ 



2H-1 



— e 2*2 



7r(s2 +a;2)2 ^/2^^ 

and this permits us to write 



2H-1 







dt x^ttV2^ dx^^ 
which concludes the proof. 

Remark 4.2 The process Iscit) = B{\C{t)\), t > has probability density 



e 2s 



^/2^7^(.s2 +<2) 



ds, X eR, t > 0. 



(4.10) 



The governing equation (jl.lip has been obtained in Nane |12 [ . 
We give here the following proof 



9^ 



t 



t 



ds 



e 2s 92 



ds 



d 



d 



t 



2I 



9s V ^2^ / as V7r(s2 +t2) 
t 



e 2s 



i9s \ V2^ j 7r(s2 +t2) 
1 92 



s=0 



ds 





e 2s 



/2^ / 7r(s2 4-^2) 



nt dx"^ 



Six) 



1 a^p 
22 a^' 



Remark 4.3 The transition law of the process {\C'^ {t)\) can be written down 
in an alternative form by exploiting the subordination relationship 



t 



7r(x2 +t2) 
Indeed we have that 



e 2s e 2s 

: t — ds 



I2^s y/2^^ 



t 



s2 + a;2 t2 + s2 



ds 



= 2 



= 2t 



e 2m e 211 

dw — 



\'2ww \/2ttw^ 

e 211, e 2z 



6 2^ e 2i 
dz — i=t- 



'2nz V2t 



/•oo 



1 



1 



\J2'KW \J2-Kz^ \j2-nz \/2wufi Jq 



ds 
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2t f°° , f°° , e~f^ e-sl 1 1 1 
dz / aw 



1 I 1 



t r r°° , 1 

= — 7 / dz / dw e 2 2 

27r^ Jo Jo w + z 

= dz [ dw „ ^ ,, e-^-^ ^ ^£;| „^ 1- (4.11) 

The result (|4.1ip shows that the density of the iterated Cauchy process still 
preserves the structure of the Cauchy distribution with exponentially distributed 
weights. 

Remark 4.4 We can give another alternative and significant representation of 

the distribution of the iterated Cauchy process. 

Since 

C{V2t)'^-—^ t>0 

we can write that 

CWC^{t)\)'= i ^, t>o. 



Furthermore by combining the previous results we can say that 

^''l^^">l''^ cMlc'(»l) - 

The last result can be also derived by observing that, for w > 

1 



> 

w 



oo roo 1_ 
* t 



A„7o 7r(.s2+a;2)7r(^ + s2) 



dsdx 



which leads to the density of , J^; ^ v a given by 

~ 7T^ W^ Jo S^ + -^^ ^+S^ TT^ Jo S^W^ + ll + tH^ 

°° t dz 2 n z t , 

zdz 



TT^ 7o Z^ + 1 1 + W TT^ Z^ + Iw'^ + t'^z'^ 

A n f 1 .^^^ ^ 2l r * y iy 



Similar calculations hold for w < 0. 
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Process 


Governing 
equation 


Iterated 
process 


Governing 
equation 


C{t) 


dq dq 
dt ~ d\x\ 

d^ 

9*2 — dx^ 


cH\cHm 


d'^q 2 I d'^q 
dt^ ~ TT^tx^ a£!;2 


BH{t) 


dp jT-L2H-ld'p 
dt — dx^ 




dq__2Hi^"-^ j^£2H-ld\ 
dt itx'^^j2Ti dx^ 


Bit) 


dp 1 d'^p 
dt ~ 2 dx^ 


Bi\cm 


d'^q 1 d^q 1 d''6 

9*2 ~ 4 dx^ TTt dx2 



Table 2 
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